Bacterial spores respond to humidity similarly to hydrogels.

Bacterial spores have outstanding properties from the materials science perspective, which allow them to survive extreme environmental conditions. Recent work by [S. G. Harrellson et al., Nature 619, 500-505 (2023)] studied the mechanical properties of Bacillus subtilis spores and the evolution of these properties with the change of humidity. The experimental measurements were interpreted assuming that the spores behave as water-filled porous solids, subjected to hydration forces. Here, we revisit their experimental data using literature data on vapor sorption on spores and ideas from polymer physics. We demonstrate that upon the change of humidity, the spores behave like rubber with respect to their swelling, elasticity, and relaxation times. This picture is consistent with the knowledge of the materials comprising the bacterial cell walls-cross-linked peptidoglycan. Our results provide an interpretation of the mechanics of bacterial spores and can help in developing synthetic materials mimicking the mechanical properties of the spores.

Ultrasonic study of water adsorbed in nanoporous glasses.

Thermodynamic properties of fluids confined in nanopores differ from those observed in the bulk. To investigate the effect of nanoconfinement on water compressibility, we perform water sorption experiments on two nanoporous glass samples while concomitantly measuring the speed of longitudinal and shear ultrasonic waves in these samples. These measurements yield the longitudinal and shear moduli of the water-laden nanoporous glass as a function of relative humidity that we utilize in the Gassmann theory to infer the bulk modulus of the confined water. This analysis shows that the bulk modulus (inverse of compressibility) of confined water is noticeably higher than that of the bulk water at the same temperature. Moreover, the modulus exhibits a linear dependence on the Laplace pressure. The results for water, which is a polar fluid, agree with previous experimental and numerical data reported for nonpolar fluids. This similarity suggests that irrespective of intermolecular forces, confined fluids are stiffer than bulk fluids. Accounting for fluid stiffening in nanopores may be important for accurate interpretation of wave propagation measurements in fluid-filled nanoporous media, including in petrophysics, catalysis, and other applications, such as in porous materials characterization.

Adsorption-Induced Deformation of Zeolites 4A and 13X: Experimental and Molecular Simulation Study.

Gas adsorption in zeolites leads to adsorption-induced deformation, which can significantly affect the adsorption and diffusive properties of the system. In this study, we conducted both experimental investigations and molecular simulations to understand the deformation of zeolites 13X and 4A during carbon dioxide adsorption at 273 K. To measure the sample's adsorption isotherm and strain simultaneously, we used a commercial sorption instrument with a custom-made sample holder equipped with a dilatometer. Our experimental data showed that while the zeolites 13X and 4A exhibited similar adsorption isotherms, their strain isotherms differed significantly. To gain more insight into the adsorption process and adsorption-induced deformation of these zeolites, we employed coupled Monte Carlo and molecular dynamics simulations with atomistically detailed models of the frameworks. Our modeling results were consistent with the experimental data and helped us identify the reasons behind the different deformation behaviors of the considered structures. Our study also revealed the sensitivity of the strain isotherm of zeolites to pore size and other structural and energetic features, suggesting that measuring adsorption-induced deformation could serve as a complementary method for material characterization and provide guidelines for related technical applications.

Surface Tension of Organophosphorus Compounds: Sarin and its Surrogates.

While the production and stockpiling of organophosphorus chemical warfare agents (CWAs), such as sarin, was banned three decades ago, CWAs have remained a threat. New approaches for decontamination and destruction of CWAs require detailed knowledge of their various physicochemical properties. In particular, surface tension is needed to describe the formation and evolution of hazardous aerosols when CWA liquids are dispersed in the air. Due to the extreme toxicity of sarin, most experimental studies are carried out using its surrogates─organophosphorus compounds which, while having similar structures, are much less toxic, e.g., dimethyl methylphosphonate (DMMP) and diisopropyl methylphosphonate (DIMP). However, not only for sarin, but also for its surrogates, literature data on the surface tension are scarce. Here we present experimental measurements and computational predictions of the surface tension of DMMP and DIMP. Classical molecular dynamics simulations using the Transferable Potentials for Phase Equilibria (TraPPE) force field produced an excellent agreement with the experimental results in the temperature range from 3 to 60 °C, validating the predictive capability of TraPPE. Consequently, we applied the TraPPE force field to sarin. Our modeled values for the sarin surface tension cover the range of temperatures from 0 to 85 °C, and the four experimental data points from the literature measured between 20 and 35 °C agree perfectly with our predictions. The temperature-dependent surface tension values for sarin and its surrogates obtained in our study can be used in models predicting the formation and evolution of aerosols made of these chemicals. Furthermore, our results justify the use of the TraPPE force field to derive the thermodynamic properties of other organophosphorus compounds with structures similar to the ones studied here.

Closer Look at Adsorption of Sarin and Simulants on Metal-Organic Frameworks.

The development of effective protection against exposure to chemical warfare agents (CWAs), such as sarin, relies on studies of its adsorption on the capturing materials and seeking candidates capable of adsorbing large amounts of sarin gas. Many metal-organic frameworks (MOFs) are promising materials for the effective capture and degradation of sarin and simulant substances. Among the simulants capable of mimicking thermodynamic properties of the agent, not all of them have been investigated on the ability to act similarly in the adsorption process, in particular, whether the agent and a simulant have similar mechanisms of binding to the MOF surface. Molecular simulation studies not only provide a safe way to investigate the aforementioned processes but can also help reveal the mechanisms of interactions between the adsorbents and the adsorbing compounds at the molecular level. We performed Monte Carlo simulations of the adsorption of sarin and three simulants, dimethyl methylphosphonate (DMMP), diisopropyl methylphosphonate (DIMP), and diisopropyl fluorophosphate (DIFP), on selected MOFs that have previously shown strong capabilities to adsorb sarin. On the basis of the calculated adsorption isotherms, enthalpy of adsorption, and radial distribution functions, we revealed common mechanisms among the particularly efficient adsorbents as well as the ability of simulants to mimic them. The findings can help in selecting a suitable simulant compound to study CWA adsorption on MOFs and guide further synthesis of efficient MOFs for the capture of organophosphorus compounds.

Molecular Simulations of Vapor-Liquid Equilibrium of Isocyanates.

The wide range of applications of the isocyanates across multiple industries sparks the interest in the study of their phase behavior. A molecular simulation is a powerful tool that can go beyond experimental investigations relying on a molecular structure of a chemical. The success of a molecular simulation relies on a description of the system, namely, force field, and its parameterization on reproducing properties of interest. In this work, we propose a united-atom force field based on the transferable potentials for phase equilibria (TraPPE) to model the vapor-liquid phase behavior of isocyanates. With Monte Carlo and molecular dynamics simulation methods and the introduced force field, we modeled vapor-liquid equilibrium for a family of linear mono-isocyanates, from methyl isocyanate to hexyl isocyanate, and hexamethylene diisocyanate. Additionally, we performed similar calculations for methyl, ethyl, and butyl isocyanates based on the all-atom GAFF-IC force field available in the literature for modeling isocyanate viscosities. We showed that the developed TraPPE-based force field generally overperformed the GAFF-IC force field and overall showed excellent performance in modeling phase behavior of isocyanates. Based on the simulated vapor pressures for the considered compounds, we estimated the Antoine equation parameters to calculate the vapor pressure in a range of temperatures. The predictions are of particular use in the investigation of thermodynamic properties for those isocyanates lacking experimental vapor pressure data. Results can also be employed in modeling the phase behavior of isocyanate mixtures to investigate their sensing and capturing. Furthermore, from the vapor-liquid equilibrium binodals, we predicted the critical properties of isocyanates which can be used in thermodynamic models based on an equation of state.

Force Fields for Molecular Modeling of Sarin and its Simulants: DMMP and DIMP.

Even three decades after signing the Chemical Weapons Convention, organophosphorus chemical warfare agents (CWAs), such as sarin, remain a threat. The development of novel methods for the detection of CWAs, protection from CWAs, and CWA decontamination motivates research on their physicochemical properties. Due to the extreme toxicity of sarin, most of the experimental studies are carried out using less toxic simulant compounds. In addition to experimental studies of sarin simulants, both sarin and simulants can be studied using in silico experiments-molecular simulations. The results of classical molecular modeling of the compounds and their agreement with experimental data rely on the force field used to describe the system. In recent years, there have been several force fields proposed for sarin and its most common simulant dimethyl methylphosphonate (DMMP). However, other simulants frequently used in experiments received less attention from the molecular simulation perspective, for example, to date, there is no force field and no simulation data for diisopropyl methylphosphonate (DIMP). Here, we compare the literature force fields for sarin and DMMP, focusing specifically on the vapor-liquid equilibrium for the pure species. We carried out Monte Carlo and molecular dynamics simulations using the existing literature force fields from which we predicted the liquid densities and vapor pressures developing the entire binodal curves. We compared the predictions to the experimental data and showed that the TraPPE-UA force field outperformed the other force fields. Thus, we extended TraPPE-UA for DIMP, utilized this force field in molecular simulations, and predicted the liquid densities and vapor pressures for a range of temperatures (binodal curve), which agreed well with the published experimental data. From the binodal, we calculated the critical properties of DIMP and demonstrated that these parameters can be used in the Peng-Robinson equation of state for this compound.

The effect of interconnections on gas adsorption in materials with spherical mesopores: A Monte Carlo simulation study.

Gas adsorption is a standard method for measuring pore-size distributions of nanoporous materials. This method is often based on assuming the pores as separate entities of a certain simple shape: slit-like, cylindrical, or spherical. Here, we study the effect of interconnections on gas adsorption in materials with spherical pores, such as three-dimensionally ordered mesoporous (3DOm) carbons. We consider interconnected systems with two, four, and six windows of various sizes. We propose a simple method based on the integration of solid-fluid interactions to take into account these windows. We used Monte Carlo simulations to model argon adsorption at the normal boiling point and obtained adsorption isotherms for the range of systems. For a system with two windows, we obtained a remarkably smooth transition from the spherical to cylindrical isotherm. Depending on the size and number of windows, our system resembles both spherical and cylindrical pores. These windows can drastically shift the point of capillary condensation and result in pore-size distributions that are very different from the ones based on a spherical pore model. Our results can be further used for modeling fluids in a system of interconnected pores using Monte Carlo and density functional theory methods.

Pore-Size Distribution of Silica Colloidal Crystals from Nitrogen Adsorption Isotherms.

Silica colloidal crystals are face-centered cubic structures comprised of silica spheres with the diameters ranging between tens and hundreds of nanometers. The voids between the spheres form pores, which can be probed by nitrogen adsorption porosimetry. Here, we prepared two mesoporous samples and a macroporous reference sample and then measured nitrogen adsorption and desorption isotherms for further characterization. We proposed a straightforward procedure for calculation of the pore-size distribution of silica colloidal crystals from nitrogen adsorption isotherms. The procedure is based on the adsorption integral equation solution with a kernel of theoretical isotherms, consistent with the procedure used for many other porous materials. The solution is carried out using the non-negative least squares (NNLS) regression with Tikhonov regularization. The kernel of mesoporous isotherms is built on the basis of the macroscopic Derjaguin-Broekhoff-de Boer (DBdB) theory of capillary condensation considering the voids as a network of spheres. Application of our procedure for the analysis of the adsorption branches of experimental isotherms resulted in bimodal distributions, where the modes matched well with the sizes of the voids in the colloidal crystals face centered cubic structure: the main mode corresponds to the octahedral voids and the second mode to the tetrahedral voids. Furthermore, we modified the surface of the samples with organics and repeated the characterization procedure for the modified samples. The resulting pore-size distribution for the samples with the modified surface matched the original one quite closely. It demonstrates the procedure as a simple and efficient technique to estimate the pore-size distribution and justifies the spherical shape approximation for the voids in the silica colloidal crystals.

In Situ Small-Angle Neutron Scattering Investigation of Adsorption-Induced Deformation in Silica with Hierarchical Porosity.

Adsorption-induced deformation of a series of silica samples with hierarchical porosity has been studied by in situ small-angle neutron scattering (SANS) and in situ dilatometry. Monolithic samples consisted of a disordered macroporous network of struts formed by a 2D lattice of hexagonally ordered cylindrical mesopores and disordered micropores within the mesopore walls. Strain isotherms were obtained at the mesopore level by analyzing the shift of the Bragg reflections from the ordered mesopore lattice in SANS data. Thus, SANS essentially measured the radial strain of the cylindrical mesopores including the volume changes of the mesopore walls due to micropore deformation. A H2O/D2O adsorbate with net zero coherent neutron scattering length density was employed in order to avoid apparent strain effects due to intensity changes during pore filling. In contrast to SANS, the strain isotherms obtained from in situ dilatometry result from a combination of axial and radial mesopore deformation together with micropore deformation. Strain data were quantitatively analyzed with a theoretical model for micro-/mesopore deformation by combining information from nitrogen and water adsorption isotherms to estimate the water-silica interaction. It was shown that in situ SANS provides complementary information to dilatometry and allows for a quantitative estimate of the elastic properties of the mesopore walls from water adsorption.

Porous Structure of Silica Colloidal Crystals.

We prepared silica colloidal crystals with different pore sizes using isothermal heating evaporation-induced self-assembly in quantities suitable for nitrogen porosimetry and studied their porous structure. We observed pores of two types in agreement with the description of silica colloidal crystals as face-centered cubic packed structures containing octahedral and tetrahedral voids. We calculated the sizes of these pores using the Derjaguin-Broekhoff-de Boer theory of capillary condensation for spherical pores. We also described the pore geometry mathematically and showed that the octahedral pore radii measured experimentally matches closely the radii of the spheres of the same volume. In the case of the tetrahedral pores, the proposed approach underestimated the pore radius by ca. 40%. Overall, this simple geometrical description provides a good representation of the porous system in silica colloidal crystals.

Mechanical Characterization of Hierarchical Structured Porous Silica by in Situ Dilatometry Measurements during Gas Adsorption.

Mechanical properties of hierarchically structured nanoporous materials are determined by the solid phase stiffness and the pore network morphology. We analyze the mechanical stiffness of hierarchically structured silica monoliths synthesized via a sol-gel process, which possess a macroporous scaffold built of interconnected struts with hexagonally ordered cylindrical mesopores. We consider samples with and without microporosity within the mesopore walls and analyze them on the macroscopic level as well as on the microscopic level of the mesopores. Untreated as-prepared samples still containing some organic components and the respective calcined and sintered counterparts of varying microporosity are investigated. To determine Young's moduli on the level of the macroscopic monoliths, we apply ultrasonic run time measurements, while Young's moduli of the mesopore walls are obtained by analysis of the in situ strain isotherms during N2 adsorption at 77 K. For the latter, we extended our previously reported theoretical approach for this type of materials by incorporating the micropore effects, which are clearly not negligible in the calcined and most of the sintered samples. The comparison of the macro- and microscopic Young's moduli reveals that both properties follow essentially the same trends, that is, calcination and sintering increase the mechanical stiffness on both levels. Consequently, stiffening of the monolithic samples can be primarily attributed to stiffening of the backbone material which is consistent with the fact that the morphology on the mesopore level is mainly preserved with the post-treatments applied.

Single Parameter for Predicting the Morphology of Atmospheric Black Carbon.

Black carbon (BC) from fuel combustion is an effective light absorber that contributes significantly to direct climate forcing. The forcing is altered when BC combines with other substances, which modify its mixing state and morphology, making the evaluation of its atmospheric lifetime and climate impact a challenge. To elucidate the associated mechanisms, we exposed BC aerosol to supersaturated vapors of different chemicals to form thin coatings and measured the coating mass required to induce the restructuring of BC aggregates. We found that studied chemicals fall into two distinct groups based on a single dimensionless parameter, χ, which depends on the diameter of BC monomer spheres and the coating material properties, including vapor supersaturation, molar volume, and surface tension. We show that when χ is small (low-volatility chemicals), the highly supersaturated vapor condenses uniformly over aggregates, including convex monomers and concave junctions in between monomers, but when χ is large (intermediate-volatility chemicals), junctions become preferred. The aggregates undergo prompt restructuring when condensation in the junctions dominates over condensation on monomer spheres. For a given monomer diameter, the coating distribution is mostly controlled by vapor supersaturation. The χ factor can be incorporated straightforwardly into atmospheric models to improve simulations of BC aging.

Molecular Simulations Shed Light on Potential Uses of Ultrasound in Nitrogen Adsorption Experiments.

Nitrogen adsorption is one of the main characterization techniques for nanoporous materials. The experimental adsorption isotherm provides information about the surface area and pore size distribution (PSD) for a sample. In this work we show that additional insight into PSD can be gained when the speed of sound propagation through a sample is measured during nitrogen adsorption experiment. We analyzed published experimental data on ultrasound propagation through a nanoporous Vycor glass sample during nitrogen adsorption experiment. Next, we calculated the change of the longitudinal and shear moduli of the sample as a function of relative vapor pressure. From this, we show that the shear modulus of the sample does not change upon filling the pores, evidencing that adsorbed nitrogen at 77 K has zero shear modulus, similarly to a bulk liquid. The longitudinal modulus of the sample behaves differently: it changes abruptly at the capillary condensation and keeps gradually increasing thereafter. We performed Monte Carlo molecular simulations to predict the compressibility of adsorbed nitrogen and then calculated the longitudinal modulus of the nitrogen-saturated Vycor using the Gassmann equation. Our theoretical predictions nicely matched the longitudinal modulus derived from the experimental data. Additionally, we performed molecular simulations to model nitrogen adsorbed in silica pores of sizes ranging from 2 to 8 nm. We found that the isothermal elastic modulus of adsorbed nitrogen depends linearly on the inverse pore size. This dependence, along with the proposed recipe for probing the modulus of adsorbed nitrogen, sets up the grounds for extracting additional information about the porous samples, when the nitrogen adsorption is combined with ultrasonic experiments.

Effect of pore geometry on the compressibility of a confined simple fluid.

Fluids confined in nanopores exhibit properties different from the properties of the same fluids in bulk; among these properties is the isothermal compressibility or elastic modulus. The modulus of a fluid in nanopores can be extracted from ultrasonic experiments or calculated from molecular simulations. Using Monte Carlo simulations in the grand canonical ensemble, we calculated the modulus for liquid argon at its normal boiling point (87.3 K) adsorbed in model silica pores of two different morphologies and various sizes. For spherical pores, for all the pore sizes (diameters) exceeding 2 nm, we obtained a logarithmic dependence of fluid modulus on the vapor pressure. Calculation of the modulus at saturation showed that the modulus of the fluid in spherical pores is a linear function of the reciprocal pore size. The calculation of the modulus of the fluid in cylindrical pores appeared too scattered to make quantitative conclusions. We performed additional simulations at higher temperature (119.6 K), at which Monte Carlo insertions and removals become more efficient. The results of the simulations at higher temperature confirmed both regularities for cylindrical pores and showed quantitative difference between the fluid moduli in pores of different geometries. Both of the observed regularities for the modulus stem from the Tait-Murnaghan equation applied to the confined fluid. Our results, along with the development of the effective medium theories for nanoporous media, set the groundwork for analysis of the experimentally measured elastic properties of fluid-saturated nanoporous materials.

Adsorption-Induced Deformation of Hierarchically Structured Mesoporous Silica-Effect of Pore-Level Anisotropy.

The goal of this work is to understand adsorption-induced deformation of hierarchically structured porous silica exhibiting well-defined cylindrical mesopores. For this purpose, we performed an in situ dilatometry measurement on a calcined and sintered monolithic silica sample during the adsorption of N2 at 77 K. To analyze the experimental data, we extended the adsorption stress model to account for the anisotropy of cylindrical mesopores, i.e., we explicitly derived the adsorption stress tensor components in the axial and radial direction of the pore. For quantitative predictions of stresses and strains, we applied the theoretical framework of Derjaguin, Broekhoff, and de Boer for adsorption in mesopores and two mechanical models of silica rods with axially aligned pore channels: an idealized cylindrical tube model, which can be described analytically, and an ordered hexagonal array of cylindrical mesopores, whose mechanical response to adsorption stress was evaluated by 3D finite element calculations. The adsorption-induced strains predicted by both mechanical models are in good quantitative agreement making the cylindrical tube the preferable model for adsorption-induced strains due to its simple analytical nature. The theoretical results are compared with the in situ dilatometry data on a hierarchically structured silica monolith composed by a network of mesoporous struts of MCM-41 type morphology. Analyzing the experimental adsorption and strain data with the proposed theoretical framework, we find the adsorption-induced deformation of the monolithic sample being reasonably described by a superposition of axial and radial strains calculated on the mesopore level. The structural and mechanical parameters obtained from the model are in good agreement with expectations from independent measurements and literature, respectively.

Modulus-pressure equation for confined fluids.

Ultrasonic experiments allow one to measure the elastic modulus of bulk solid or fluid samples. Recently such experiments have been carried out on fluid-saturated nanoporous glass to probe the modulus of a confined fluid. In our previous work [G. Y. Gor et al., J. Chem. Phys., 143, 194506 (2015)], using Monte Carlo simulations we showed that the elastic modulus K of a fluid confined in a mesopore is a function of the pore size. Here we focus on the modulus-pressure dependence K(P), which is linear for bulk materials, a relation known as the Tait-Murnaghan equation. Using transition-matrix Monte Carlo simulations we calculated the elastic modulus of bulk argon as a function of pressure and argon confined in silica mesopores as a function of Laplace pressure. Our calculations show that while the elastic modulus is strongly affected by confinement and temperature, the slope of the modulus versus pressure is not. Moreover, the calculated slope is in a good agreement with the reference data for bulk argon and experimental data for confined argon derived from ultrasonic experiments. We propose to use the value of the slope of K(P) to estimate the elastic moduli of an unknown porous medium.

Deformation of Microporous Carbons during N2, Ar, and CO2 Adsorption: Insight from the Density Functional Theory.

Using the nonlocal density functional theory, we investigate adsorption of N2 (77 K), Ar (77 K), and CO2 (273 K) and respective adsorption-induced deformation of microporous carbons. We show that the smallest micropores comparable in size and even smaller than the nominal molecular diameter of the adsorbate contribute significantly to the development of the adsorption stress. While pores of approximately the nominal adsorbate diameter exhibit no adsorption stress regardless of their filling level, the smaller pores cause expansive adsorption stresses up to almost 4 GPa. Accounting for this effect, we determined the pore-size distribution of a synthetic microporous carbon by simultaneously fitting its experimental CO2 adsorption isotherm (273 K) and corresponding adsorption-induced strain measured by in situ dilatometry. Based on the pore-size distribution and the elastic modulus fitted from CO2 data, we predicted the sample's strain isotherms during N2 and Ar adsorption (77 K), which were found to be in reasonable agreement with respective experimental data. The comparison of calculations and experimental results suggests that adsorption-induced deformation caused by micropores is not limited to the low relative pressures typically associated with the micropore filling, but is effective over the whole relative pressure range up to saturation pressure.

Adsorption-Induced Surface Stresses of the Water/Quartz Interface: Ab Initio Molecular Dynamics Study.

Adsorption-induced deformation is expansion or contraction of a solid due to adsorption on its surface. This phenomenon is important for a wide range of applications, from chemomechanical sensors to methane recovery from geological formations. The strain of the solid is driven by the change of the surface stress due to adsorption. Using ab initio molecular dynamics, we calculate the surface stresses for the dry α-quartz surfaces, and investigate how these stresses change when the surfaces are exposed to water. We find that the nonhydroxylated surface shows small and approximately isotropic changes in stress, while the hydroxylated surface, which interacts more strongly with the polar water molecules, shows larger and qualitatively anisotropic (opposite sign in xx and yy) surface stress changes. All of these changes are several times larger than the surface tension of water itself. The anisotropy and possibility of positive surface stress change can explain experimentally observed surface area contraction due to adsorption.